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Abstract 

The six quark state(uuddss) called H dibaryon(J p = + ,S = —2) has 
been calculated to study its existence and stability. The simulations are 
performed in quenched QCD on 8 3 x 24 and 16 3 x 48 anisotropic lattices 
with Symanzik improved gauge action and Clover fermion action. The 
gauge coupling is /3 = 2.0 and aspect ratio £ = a s /at = 3.0. Preliminary 
results indicate that mass of H dibaryon is 2134(100)Mev on 8 3 x 24 lattice 
and 2167(59)Mev on 16 3 x 48 respectively. It seems that the radius of H 
dibaryon is very large and the finite size effect is very obvious. 

1 Introduction 

In 1976, Jaffc pointed out that the quark bag model predicted the existence of 
H dibaryon, which is a compound state of 6 quarks (uuddss) [I]. It is the lowest 
bound state in dibaryon sector and will be a spin strangeness -2 SU(3) flavor 
singlet. Jaffe's original bag model suggested that tuh = 2150Mev, 81Mev below 
the 2231Mev AA threshold. 

Since Jaffe's prediction, people tried to find out the H-Dibaryon state both 
by experiment and theoretical calculation. On the experimental side, until now, 
there are no enough evidences on the existence of the H-Dibaryon. On the 
theoretical side, many theoretical calculations have been maked to predict the 
mass of H dibaryon. One of the most efficient ways to study this state is from the 
first principle of QCD, i.e Lattice QCD. Many simulation results suggested that 
H dibaryon is a bound state, but some gave contrary conclusions [2 O 01 [6] . 

We perform the numerical simulations with refined methods, which includes 
improved gauge and fermion action, smearing techniques and tadpole improve- 
ment. To reduce computer cost and determine large masses more accurately, 
we do the simulations on anisotropic lattice for both gauge action and fermion 
action. This is our advantage comparing with what the forthgoer have done. 
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2 H-Dibaryon Simulation Details 

2.1 Actions 

We generate the configurations using improved, anisotropic action[7]: 

q n r ^ ^sp . 4 £^£p 1 ^sr 1 C^str -i \ 

^ ~ W 3R + " 12 R ~ 12^p UJ 

where (3 = 6/<? 2 , g is the QCD coupling, u s and Ut are the mean-link renor- 
malization parameters, £ is the aspect ratio (£ = a s /at at the tree level in 
perturbation theory), and 



a *P = EE ^Tr[l - ^(1)^(1 + i)C/?(a; + j)U}(x)], 

x i>j 

= E E \ ReTr \ l - Ui{x)U t {x + i)U}{x + t)U}(x)}, 

x i 

n «" = EE ^eTr[l - Ui(x)Ui(x + i)U s (x + 2i)Zl](x + j + i)U}(x + j)Uj(x)}, 



X ijij 



n-r = E E >7Mi - + w . + * )t /? ( , + * + i)u} ( X + t) ^( a) ], 

where x labels the sites of the lattice, i,j are spatial indices, and U^{x) is the 
parallel transport matrix in the gauge field from site x to x + /i. 

For the quark action, we employ the space-time asymmetric clover quark 
action on anisotropic lattice 8 9] [10 \ 

s f = e* a 

E Et*^( rs _ 7i)Ci,x*«-H + **(r s + 7i)tfjx-<*x-i] 

-^ t x;[^( i -7t)^,»*»+t+*»(i+7t)^Vt**-*] 

+iiT s c s E + iif s ct E ^xO-ti-PiiCa:)^, (2) 

x,i<j x,i 

where K Sit and c Sj t are the spatial and temporal hopping parameters and the 
clover coefficients, respectively. The hopping parameters K s>t are related to the 
bare quark mass mo = atm q o through 

a t m q0 ee l/{2K t ) - 3r s /( - 1, C = (3) 

We perform the simulations using tree-level improved Symanzik action and 
Clover fermion action, with gauge coupling j3 = 2.0 and the aspect ratio £ = 3.0. 
The lattice sizes are 8 3 x 24 and 16 3 x 48. 
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2.2 Operator and Correlation Function of H Dibaryon 

The operator of H-Dibaryondl : 

Oh(x) = 3(udsuds) — 3(ussudd) — 3(dssduu), (4) 

(abcdef) = e abc e de f{Cj 5 ) a/3 {Cj 5 ) lS (Cj 5 ) e4> a*(x)b b p(x)c c t (x)d*(x)e e s (x)fl{x), 

(5) 

where a, b, c, d, e, / are color indices and a, (3, 7, S, e, (f> are spinor indices. 

And the corresponding correlation function of H Dibaryon can be written as 
Gh{x,t) =< Oh(x, t)Ojj(O) >, which involves terms of the structure [T2]: 

(UnU 22 - UuUnXDnDm - D 12 D 21 )(S n S 2 2 - S X2 S 2X ). (6) 

To decide whether the H Dibaryon is stable or not, usually we can compare 
the mass of AA with H-Dibaryon's. The AA operator [12]: 

A (x) = eabc (C 75 )^[<(x)4(x)4( : r) + <(a;)4(x)^( : r)-2<(ar)4(ar)^( a; )]. 

(7) 

2.3 Smearing Techniques 

To reduce the excited-state contamination, we use the smearing techniques 
which can provide a better overlap with the ground state [13] ■ For the quark 
fields we use: 

ip'(x,R) = ^2(U\x-fi)...U ji {x-Rfi)Tp{x-Rfi) 

+U(x) . . . U(x + (R- l)p)i){x + RfM)). (8) 

A more large plateau in the region with small errors is obtained with smear- 
ing. 

3 Simulation Results 

Out results are shown in Tables 1, 2, 3 and 4. The kjj t is the temporal heavy 
kappa, which corresponds to s quark; the fc^ t is temporal light kappa, which 
corresponds to u and d quarks. 

The H-Dibaryon remain lighter than two A at all combinations of the hopping 
parameters, as shown in tables. 

To obtain the physical masses of H and A, one has to extrapolate or interpo- 
late the fci, t and kjj t to physical hopping parameters. Since (m w a) 2 is linearly 
related to 1/k, we can determine the critical hopping parameter k c at which 
(m^a) 2 vanishes. We take physical k u d as k c because they are very close. The 
physical k s can be determined from the ratio of a strangeness carrying particle 
to a non-strange one, here we obtain the k s by the mass ratio of lambda and 
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Table 1: ma t of the A (8 3 x 24,(3 = 2.04 = 3-0) 



k Ht 1 k Lt 


0.23810 


0.23923 


0.24039 


0.24155 


0.24272 


0.23256 


1.3772(65) 


1.3496(67) 


1.3202(69) 


1.2885(71) 


1.2542(73) 


0.23365 


1.3665(66) 


1.3388(68) 


1.3093(70) 


1.2777(72) 


1.2433(74) 


0.23474 


1.3553(66) 


1.3276(68) 


1.2981(70) 


1.2664(72) 


1.2319(75) 


0.23586 


1.3438(67) 


1.3160(69) 


1.2864(71) 


1.2546(73) 


1.2201(75) 


0.23697 


1.3317(68) 


1.3039(69) 


1.2742(71) 


1.2424(74) 


1.2078(76) 


0.23810 


1.3189(69) 


1.2910(71) 


1.2613(73) 


1.2294(75) 


1.1947(77) 



Table 2: ma t of the A (16 3 x 48,/? = 2.04 = 3-0) 



knJk Lt 


0.23810 


0.23923 


0.24039 


0.24155 


0.24272 


0.23256 


1.4161(32) 


1.3881(33) 


1.3583(34) 


1.3262(34) 


1.2913(35) 


0.23365 


1.4052(32) 


1.3772(33) 


1.3473(34) 


1.3152(34) 


1.2802(35) 


0.23474 


1.3940(33) 


1.3659(33) 


1.3359(34) 


1.3037(35) 


1.2688(36) 


0.23585 


1.3822(33) 


1.3541(34) 


1.3241(34) 


1.2918(35) 


1.2568(36) 


0.23697 


1.3699(33) 


1.3418(34) 


1.3117(34) 


1.2794(35) 


1.2443(36) 


0.23810 


1.3568(32) 


1.3286(33) 


1.2985(33) 


1.2661(34) 


1.2309(35) 



nucleon. Our calculations suggest that k c — 0.25256(41) and k s — 0.2413(40) 
on 8 3 x 24 lattice and k c = 0.25323(18) and k s = 0.2422(21) on 16 3 x 48 lattice. 

In fig.l, we performed a linear fit to extrapolate the ran and m\ to the 
physical k s , and obtained the H's mass tuh = 2134(100)Mev, which is lower than 
two A's. The difference in mass is tuh ~ 2m\ — — 97(100)Mev. That means the 
H-Dibaryon tends to be a bound state but actually we can't make this conclusion 
because the error is larger than the mass difference on 8 3 x 24 lattice. In fig. 2, 
on larger lattice(16 3 x 48), m H = 2167(59)Mev and m H - 2m A = -64(59)Mev, 
indicates that the energey of H dibaryon tends to below the AA threshold. 
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Tabic 3: ma t of the H (8 3 x 24,/3 = 2.0,£ = 3.0) 



k Ht 1 k Lt 


0.23810 


0.23923 


0.24039 


0.24155 


0.24272 


0.23256 


2.656(17) 


2.603(18) 


2.545(18) 


2.484(19) 


2.417(19) 


0.23365 


2.635(18) 


2.581(18) 


2.523(18) 


2.462(19) 


2.395(20) 


0.23474 


2.613(18) 


2.558(18) 


2.501(19) 


2.439(19) 


2.372(20) 


0.23586 


2.589(18) 


2.535(18) 


2.477(19) 


2.415(19) 


2.348(20) 


0.23697 


2.565(18) 


2.510(19) 


2.452(19) 


2.390(20) 


2.323(20) 


0.23810 


2.539(18) 


2.484(19) 


2.426(19) 


2.364(20) 


2.297(20) 



Table 4: ma t of the H (16 3 x 48,/3 = 2.0,£ = 3.0) 





0.23810 


0.23923 


0.24039 


0.24155 


0.24272 


0.23256 


2.7641(87) 


2.7093(90) 


2.6508(93) 


2.5882(97) 


2.5203(102) 


0.23365 


2.7422(88) 


2.6872(91) 


2.6287(94) 


2.5660(98) 


2.4980(103) 


0.23474 


2.7194(89) 


2.6643(92) 


2.6057(95) 


2.5429(99) 


2.47489(105) 


0.23585 


2.6956(90) 


2.6405(93) 


2.5818(97) 


2.5188(101) 


2.4507(106) 


0.23697 


2.6708(91) 


2.6155(94) 


2.5567(98) 


2.4937(102) 


2.4254(108) 


0.23810 


2.6407(86) 


2.5854(88) 


2.5265(92) 


2.4634(96) 


2.3951(101) 




II 

AA 



4.45 



13=2.0 i;=3.0 V0l=8 3 x24 
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p=2.0 E=3.0 vol=16x48 



4 Conclusion and Future Plans 

We have presented the results of the lattice investigation on the H dibaryon state 
employing anisotropic improved gauge and anisotropic clover fcrmion actions. 
The advantage of using anisotropic QCD actions is to get get a better signal so 
as to obtain a large plateau on small lattice size which can save us much more 
computer cost. In the mean time, the simulation results are more accurately. 

Our results indicate that, both on the 8 3 x 24 and 16 3 x 24 anisotropic 
lattice, the masses of H dibaryon are less than that of two As. It seems that the 
H dibaryon do exist as a bound state. 

The masses of H Dibaryon on two different lattices are not so close and we 
believe that the finite size effect of H dibaryon should be taken into account. 
We plan to calculate the H dibaryon on larger lattice size to further study the 
finite size effect. We also intend to calculate other six quarks states, such as 
possible proton-antiproton, deuteron, and so on. It should be very interesting. 
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